SKYRME-RPA DESCRIPTION OF SPIN-FLIP Ml GIANT RESONANCE 
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The spin-flip Ml giant resonance is explored in the framework of Random Phase Approximation on 
the basis of the Skyrme energy functional. A representative set of eight Skyrme parameterizations 
(SkT6, SkM*, SLy6, SG2, SkO, SkO', SkI4, and SV-bas) is used. Light and heavy, spherical 
and deformed nuclei ( 48 Ca, 158 Gd, 208 Pb, and 238 U) are considered. The calculations show that 
spin densities play a crucial role in forming the collective shift in the spectrum. The interplay of 
the collective shift and spin-orbit splitting determines the quality of the description. None of the 
considered Skyrme parameterizations is able to describe simultaneously the Ml strength distribution 
in closed-shell and open-shell nuclei. It is found that the problem lies in the relative positions of 
proton and neutron spin-orbit splitting. Necessity to involve the tensor and isovector spin-orbit 
interaction is called for. 



I. INTRODUCTION 

Nuclear density- functional-theory (DFT), with the 
most prominent representatives being Skyrme-Hartree- 
Fock (SHF), the Gogny forces, and the relativistic mean- 
field model, achieved a high level of quality in the de- 
scription of ground state and dynamics of atomic nuclei 
[E 0, 9 ■ Most applications for nuclear dynamics up to 
now have been concerned with electrical excitation modes 
(natural parity). The description works generally well, 
except for some persistent problems with the isovector gi- 
ant resonance (GR) in light nuclei [HQ. Much less work 
has been done yet for magnetic excitations (unnatural 
parity). At the same time, magnetic modes are sensi- 
tive to a different class of force parameters, namely those 
related to spin. An exploration of magnetic resonances, 
like spin-flip Ml, could essentially improve the spin-orbit 
interaction. Magnetic modes could clarify the role of still 
vague tensor interaction Q. Also, the spin-flip Ml reso- 
nance is a counterpart of the spin-isospin Gamow- Teller 
resonance which is of great current interest in connection 
with astrophysical problems 0, H, 0, H, Q • Investigation 
of the Ml resonance could be useful in this connection as 
well. 

There are many studies of the spin-flip Ml mode within 
simple models beyond the DFT, see e.g. reviews [T(| 
EH, 112] • At the same time, as far as we know, a DFT 
treatment is limited to a few publications using SHF [13, 
[HI and even that is not carried through fully consistently. 
The work [l3| uses a hybrid model with partial inclusion 
of SHF in the Landau-Migdal formulation [r| ■ The other 
study [HI uses the early Skyrme forces and omits the 
crucial spin density. These studies, being a useful first 
step, are not satisfactory for nowadays demands. 

The present study aims at a fully self-consistent de- 
scription of the spin-flip Ml mode in the framework of 
SHF. Previous investigations for Gamow- Teller @, H, [1| 
and spin-flip Ml modes [l3| hint that spin-density re- 



sponse could be decisive to get a sizeable collective shift. 
So all the Skyrme terms with spin density (usually omit- 
ted in calculations for electric modes) have to be im- 
plemented and scrutinized. Furthermore, because of the 
obvious importance of spin-orbit splitting, the responses 
delivered by spin-orbit and tensor interactions have to 
be inspected as well. Since the quality of the description 
may depend on the particular Skyrme parameterization 
as well as on nuclear shape and mass region, a variety of 
parameterizations should be checked for light and heavy, 
spherical and deformed nuclei. Note that the Ml mode 
in heavy open-shell nuclei (rare-earth and actinides) ex- 
hibits a pronounced double-peak structure [HI, HU while 
closed-shell nuclei ( 48 Ca, 208 Pb, ...) show only one peak 
[l7l [Lsl [l9| . All these demands are met in the present 
study which, to the best of our knowledge, is the first sys- 
tematic and self-consistent SHF exploration of spin-flip 
Ml. 



II. DETAILS OF CALCULATIONS 



We will consider the resonance in the doubly-magic 
nuclei 48 Ca and 208 Pb and axially deformed nuclei 158 Gd 
and 238 U. A representative set of eight SHF parameteri- 
zations is used: SkT6 (H, SkO [H[, SkO' [H[, SG2 [H, 
SkM* 23], SLy6 [H, SkI4 H, and SV-bas (|. They 
exhibit a variety of effective masses (from m*/m=l in 
SkT6 down to 0.65 in SkI4) and other nuclear matter 
characteristics. Some of the forces (SLy6) were found fa- 
vorable in the description of E1(T=1) GR [2H WH [2gj [29| . 
Others were used in studies of Gamow- Teller strength 
(SG2, SkO') ©, H3] or peculiarities of spin-orbit 

splitting (SkM) j25[. The forces SkT6, SG2 and SkO' in- 
volve the tensor spin-orbit term. The parameterization 
SV-bas represents one of the latest upgrades in SHF. 

The calculations are performed within the self- 
consistent separable random-phase-approximation 
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(SRPA) ap pro ach based on factorized Skyrme residual 
interaction |26l . l30l l3l| . The self-consistent factorization 
considerably reduces the computational expense of RPA 
while maintaining a high accuracy. This allows to per- 
form systematic studies in both spherical and deformed 
(heavy and super-heavy) nuclei HE \2]l IH, US E}] • The 
residual interaction includes all contributions arising 
from the SHF functional as well as the Coulomb (direct 
and exchange) and pairing (at BCS level) terms [2g.[3l|). 
The Skyrme energy density to be exploited reads 0, Q 
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where frj, 6^, &^ are the force parameters. This func- 
tional involves time-even (nucleon p q , kinetic-energy r q , 
spin-orbit 3 q ) and time-odd (current j q , spin s q , vec- 
tor kinetic-energy T q ) densities where q denotes pro- 
tons and neutrons. Densities without index, like p = 
Pp + Pni denote total densities. The contributions with bi 
(i=0,l,2,3,4) and b[ (i=0,l,2,3) are the standard terms re- 
sponsible for ground stateproperties and electric excita- 
tions of even-even nuclei [JQ. In the standard SHF, the 
isovector spin-orbit interaction is linked to the isoscalar 
one by b' A = 64. The tensor spin-orbit terms oc b\, b\ are 
often skipped. In (TTJ) they can be switched by the param- 
eter 7t. The spin terms with bi, b\ become relevant only 
for odd nuclei and magnetic modes in even-even nuclei. 
Though bi,b'i may be uniquely determined as functions 
of bi,b[ 3], their values were not yet well tested by nu- 
clear data. Moreover, following a strict DFT, they can 
be considered as free parameters. Just these spin terms 
may be of a paramount importance for the spin-slip Ml. 
Hence all them are taken into account in SRPA. 

In addition to second functional derivatives entering 
the SRPA residual interaction for electric modes, 

S 2 E S 2 E S 2 E 5 2 E ^ 
8p q ,5p q ' 5r q ,5p q ' SJ q >5p q ' 6j q ,6j q 



the present treatment also involves the terms with 

S 2 E S 2 E 6 2 E S 2 E 
Sj q ,Ss q ' 5s q >5s q ' S3 q ,5J q ' 5T q <6s q ' 

SRPA generators include spin and orbital input opera- 
tors P° = R(r)s q + and P l q = R(r)l q + with R(r) being 1 
or r 2 . In deformed nuclei, to take into account the cou- 
pling between spin and quadrupole = 1 + states, the 
quadrupole generator Q q = r 2 Y2i(0) is also added. The 
convergence of the results with including more generators 
was checked. See details in [32| . 

The SHF calculations employ a coordinate-space grid 
with the mesh size 0.7 fm. For deformed nuclei, cylindri- 
cal coordinates are used and the equilibrium quadrupole 
deformation is found by minimization of the total energy 
[26L I29I] . The single-particle spectrum involves all lev- 
els from the bottom of the mean field well up to +20 
MeV. In the heaviest nucleus under consideration, 238 i7, 
this results in ~17000 two-quasiparticle (2qp) K* = 1+ 
pairs with the excitation energies up to 50-70 MeV. Note 
that for electric E1(T=1) and E2(T=0) excitations such 
single-particle space provides a satisfying exhaustion of 
the energy- weighted sum rules [28| . 

The spectral distribution of the spin-flip Ml mode with 
K v = 1 + is presented as the strength function 

S{Ml-u) = |(* v |M|*o)| 2 C(w - w„) (4) 

where — lj v ) — A/[27r[(w — uj v ) 2 + ^-]] is a Lorentz 
weight with the averaging parameter A=l MeV. Such av- 
eraging width is found optimal for the comparison with 
experiment and simulation of broadening effects beyond 
SRPA (escape widths, coupling with complex configura- 
tions). Further, \l/o is the ground state, v runs over the 
RPA K w = 1+ states with energies uj v and wave func- 
tions \I/„. The operator of spin- flip Ml transition reads 

in standard notation as M = P'B\f^Ylqlds^+ + 9?^+\ 
with spin g-factors g% — 5.58<j p and = — 3.82<;„ 
quenched by <r p =0.68 and ^„=0.64. As we are inter- 
ested in the spin-flip Ml, the orbital response is omitted, 
i.e. we put q q = 0. Note that in the experimental data 
[HI [H, Tl7l . ll8L [HI used for the comparison, the orbital 
contribution is strongly suppressed. The strength func- 
tion ([?]) is computed directly, i.e. without calculation of 
RPA states v, which additionally reduces the computa- 
tion effort Ha H ED, 111 • 

More details of SRPA formalism are given in the ap- 
pendices A, B, and C. 



III. RESULTS AND DISCUSSION 

In Fig. 1 the collective shifts of the main resonance 
peak in 48 Ca, 208 Pb, 158 Gd, and 238 U, obtained with dif- 
ferent Skyrme parameterizations, are shown. The shifts 
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FIG. 1: Collective shifts of the Ml peak for different Skyrme 
parameterizations (as indicated along the x-axis) in 48 Ca, 
208 Pb, 158 Gd and 238 U. The plots exhibit the total (black 
boxes) shifts as well as the partial ones with bo and 63 (open 
circles), 62 (open triangles), and 61 (stars). The 61 contribu- 
tion exists only for SkT6, SG2, and SkO'. For better view the 
symbols are connected by lines. 



are defined as E s hi/t = Esrpa — E^ qp , i.e. as a dif- 
ference in the energies of SRPA and unperturbed two- 
quasiparticle Ml peaks. The 2qp strength is calculated 
by using without the residual interaction. In addition 
to the total shift, the contributions from different spin- 
density dependent terms as well as from the tensor force 
(for SkT6, SG2 and SkO') are shown. It is seen that the 
total collective shifts are generally modest and vary from 
1-2 MeV in 48 Ca to 0.5-2 MeV in 208 Pb and 0.5-1.5 MeV 
in 158 Gd and 238 U. The low value emerges from contri- 
butions pulling in different directions. This holds for the 
separate shifts from 60 and 63 (not disentangled here). 
The &2-term gives a negative shift in contrast to the pos- 
itive one from bo, 63. Anyway the contribution from bo, &3 
usually dominates thus giving the total upshift in accor- 
dance with isovector character of the resonance. All the 
forces give generally similar results. Note a sizable con- 
tribution of the tensor interaction for SkT6 and SG2. 
For SkO' this contribution is negligible, except for 48 Ca 
where it is so strong that gives a negative total E s hift- 
It should be emphasized that the non-spin contributions 
(with 6j, ty) alone do not provide any collective shift and 
leave the Ml strength unperturbed. The whole shift is 
produced by the spin-dependent terms oc bi, b^. 

The calculations give a reasonable summed B(M1) 
strength. In the interval 0-45 MeV the unperturbed 
strength is 3.2 and 18.4-18.6 fi% in 48 Ca and 208 Pb. The 
residual interaction changes these values and we have 2.5 
- 4.8 fi 2 N in 48 Ca and 14.8 - 17.3 fi% in 208 Pb as compared 
with experimental values ~ 5.3 [i 2 N [llj and ~ 17.9 ^i 2 N 
[HI], respectively. Note a strong collective effect in 48 Ca. 

It is well known that proton and neutron spin-orbit 
splittings E1 a represent another crucial ingredient in de- 




FIG. 2: a): The spin-flip Ml resonance in 238 U calculated 
with SG2 force: total (bold solid line), proton (dotted line) 
and neutron (solid line) strengths. The experimental data 
are given by black boxes with bars. b)-d): proton and neutron 
unperturbed spin-orbit splittings for different Skyrme force in 

158 Gd 238TJ and 208p b _ 



48 



Ca E =10.2 MeV 
exp 

SG2 SV-bas I 1 ° 
SkO' \n a SLy6 




208 nu E =7.3 MeV 

rb <»j 



SkO' 

V 



SV-bas 



12 14 4 

E [MeV] 

FIG. 3: The spin-flip Ml resonance in 48 Ca (left) and 208 Pb 
(right) calculated within SRPA for 8 Skyrme forces as indi- 
cated. For better distinguishability, the strength in 208 Pb for 
SG2, SLy6, SkT6 and SkM* is depicted by the bold line. The 
experimental energies in 48 Ca ;i7l and 208 Pb IISH are marked 
by vertical arrows. 



scription of spin-flip Ml mode [lfj, lll|, |l2j. Usually 
E p so < E™ which leads in rare-earth and actinide nu- 
clei to a two-peak structure of the resonance with domi- 
nant proton (neutron) origin of the lower (upper) peak. 
This is demonstrated for 238 U in Fig. 2a) where the pro- 
ton and neutron components of the Ml mode, obtained 
with <^=0.68, <;„=0 and <? p =0, ^„=0.64, respectively, are 
shown. Panels b)-d) exhibit proton and neutron split- 
tings E] a for different Skyrme forces. The splittings are 
evaluated from centroids of the proton and neutron peaks 
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FIG. 4: The spin-flip Ml resonance in 58 Gd described with 
eight Skyrme parameterizations as indicated. The experimen- 
tal data are from [TBI, [lrj]. 



of the unperturbed Ml strength. The results strongly de- 
pend on the parameterization. In most of the cases we 
have EP so < with E% = |£& - E^\ -1-2 MeV but 
SkI4, SkO and SkO' give very close splittings with even 
EP a > E™ for SkI4 and SkO. The latter is related with a 
low value of 64 and nonzero value of b' 4 in SkI4. Note that 
the experimental proton-neutron splitting in Ml gross 
structure is ~ 2 MeV in 158 Gd and 238 U [H, 13 and 
zero in 208 Pb [H[l|. 

In Fig. 3 the SRPA Ml strength function (gj) in spher- 
ical doubly-magic nuclei 48 Ca and 208 Pb is presented. In 
48 Ca, the resonance is produced only by neutron spin- flip 
transition u(lfy 2 , 1/5/2) yielding the one-peak structure. 
This feature is correctly reproduced by all the parame- 
terizations. However, most of them underestimate the 
resonance energy (with worst SkO' case because of the 
strong and possibly wrong tensor contribution) and only 
the forces SLy6, SkI4, and SkM* (with maximal E™ ) 
give the Ml energy close to experiment. The success of 
these forces is obviously determined by a suitable neutron 
spin-orbit splitting. 

However, the same figure shows that in 208 Pb the forces 
SLy6, SkM, and SkM* considerably overestimate the Ml 
energy while the best result is achieved by SkO. Note 
that only SkO, SkO' and SkM, all having a small 
give a one-peak resonance structure in accordance with 
experiment [HI]. This is because only for these parame- 
terizations the interaction energy ( w collective shift) is 
larger than E™£ and so a significant mixture of proton 
and neutron components with forming of a one-peak res- 
onance becomes possible. The other forces have too large 
E™g and produce the two-peak structure. This demon- 
strates the great importance of the interplay between the 
residual interaction and relative proton and neutron spin- 
orbit splitting E™p for the description of spin-flip Ml. 

Figs. 4 and 5 present SRPA results for deformed 158 Gd 
and 238 U. Here, in contrast to 208 Pb, the experiment 
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FIG. 5: The same as in Fig. 4 for 2JS XJ. 



yields a double-peak structure and so the one-peak pic- 
ture from SkO, SkO', and SkM fails. The description is 
generally quite poor, with exception of SV-bas in 158 Gd 
and SG2 in 238 U. Thus we see that any Skyrme param- 
eterization fails to describe simultaneously the one-peak 
structure in closed-shell nuclei and two-peak structure in 
open-shell nuclei. The reason is yet unclear. However, 
the relative spin-orbit splitting E"P is evidently one of 
the key factor and it is important to find a way to con- 
trol it. 

Fig. 6 explores the influence of the spin-orbit contribu- 
tions by systematic variation of tensor spin-orbit (panel 
a) and isovector spin-orbit terms (panel b). Variation 
of the tensor spin-orbit strength shifts significantly both 
E" a and E P a while leaving the relative order unchanged. 
But the variation of the isovector spin-orbit strength 6 4 
has a strong effect on the relative positions. So by simul- 
taneous monitoring tensor and isovector spin-orbit inter- 
actions one may control better the spin-orbit splittings. 
Besides the strong effect on single-particle energies, these 
interactions also affect the collective shifts (see Fig. 1). 
Altogether they represent very promising tool for further 
improvement of the description of spin-flip Ml modes. 



IV. CONCLUSIONS 

We have studied the ability of Skyrme forces to de- 
scribe the spin-flip Ml resonance using RPA with self- 
consistent factorized residual interaction. The results 
show that the terms with spin and spin-orbit densities 
are responsible for a sizable collective shift of the reso- 
nance peak. The spin-orbit splitting of the underlying 
two-quasiparticle states is of crucial importance for the 
final pattern of the spectrum (single-peak versus double- 
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FIG. 6: Dependence of unperturbed spin-orbit splitting for 
proton 7r(/i~^ 2 , /i g / 2 ) and neutron v{i~^ 2 ,i 11 f 2 ) configura- 
tions in 208 Pb on (a) the attenuation < jt < 1 for tensor 
interaction and (b) the parameter b' 4 of isovector spin-orbit 
interaction. Skyrme parameterizations with varied 7t or i> 4 
are used [jj. The proton (neutron) splittings are marked by 
circles (boxes). Filled symbols mark results for parameteri- 
zations refitted for given 7t or 64. Open symbols stand for 
SV-bas where only 7t or & 4 are varied, i.e. without refitting. 



peak structure). The residual interaction tends to mix 
of proton and neutron spin-orbit partners and so works 
towards a one-peak structure, as in 208 Pb. However, a 
large difference between neutron and proton spin-orbit 
splitting inhibits this mixing and produces two distinct 
proton and neutron peaks, as seen experimentally in rare- 
earth and actinide nuclei. 

None of eight Skyrme parameterizations used in the 
present study is able to describe simultaneously the one- 
peak and two-peak structures in closed-shell and open- 
shell nuclei. In most of the cases, the resonance energies 
are badly reproduced as well. A first exploration indi- 
cates that fine-tuning of tensor and isovector spin-orbit 
interactions could improve the description, as they af- 
fect both the spin-orbit splittings and collective shifts. 
Work in this direction is in progress. A corresponding 
improvement of the Skyrme parameterizations would be 
important not only for description of spin-flip Ml res- 
onance (which is a challenge itself) but also for better 
treatment of the spin-orbit interaction in nuclei. 
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APPENDIX A: SRPA 

The SRPA formalism is given elsewhere [26|, [3(J l3l| - 
So here we will present only principle points and pecu- 
liarities pertinent to magnetic excitations. The SRPA 
approximates the residual interaction of Skyrme RPA in 
a factorized (separable) form as 

. 1 K 

Kcs P — -77 X/ X/ { K qk,q'k'X qk X q/kr +rj qk ^ k ,Y qk Y q , k >} 

qq' k 7 k' = l 

( A1 ) 

where the indices q and q' label neutrons and protons, 
k numbers the separable terms, X qk and Y qk are time- 
even and time-odd hermitian one-body operators, K qkjq ' k i 
and rj qk ^ q ' k i are the corresponding strength matrices. We 
need these two kinds of the operators since the relevant 
Skyrme functionals involve both time-even and time-odd 
densities, see 0, H, M, U HI . 

The starting point is a Skyrme functional E[ J q (r, t)] = 
J <ir7-fek(r, t), say that in ([1]), with a set of local densities 
J q . The separable operators and strength matrices are 
derived from the functional by using the scaling transfor- 
mation for the perturbed wave function of the system: 

K 

!*(*))« = Y\_^[-iQqk{t)Pqk]exp[-ip qk (t)Q qk ]\) q 
k=l 

(A2) 

were both \^f(t]) q and ground state |) g are Slater deter- 
minants, <9 g fc(r) and Pqfc(r) are generalized coordinate 
(time-even) and momentum (time-odd) hermitian one- 
body input operators, H stands for the full Hamilto- 
nian. For EA modes, the input operators Q g fe(r) deter- 
mine -Pqfc(r) = i[H, Q q k], while for MA modes, the input 
operators P g fe(r) gives Q qk (r) = i[H,P qk }. Further 

Qqk(t) = Qqk cos(ujt) ,p q k(t) = p qk sin(wt) (A3) 

are corresponding collective variables. As seen from (|A2|) . 
the input operators are determined up to arbitrary con- 
stant multipliers whose action can be compensated by 
the proper rescaling the collective variables. 

The number K of input operators in (|A2[) determines 
the number of the separable terms in V^P . The treat- 
ment converges to exact RPA for K — > 00. In practice, 
a good approximation to RPA is already obtained for a 
small K = 2 -j- 5 if the input operators Q qk (P q k) are 
properly chosen. 

The separable operators and (inverse) strength matri- 
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ces in (|A1|) are self-consistently constructed as [2a. I3CI l3ll 



a'aq' <? *? 



[P gfe) >])Jf5(A4) 



q' a'aq' 9 9 



g' k' ,qk 



»7. 



-l 

q' k' ,qk 



6 2 E 



where J q are the density operators. 
The final RPA equations have the form 

qk 

Y.^ F %^ k +P^ F %V qk -%,l. qk ]} = 0(A9) 

qk 



with 



r q'k',qk 



2 E aAS 

ij" ,ph£q" 



{ph\A* k \y(ph\&; k ,\} 



e 1h - u i 



and 



e p/l , for A = B 
ctAB = | — it^is, for A = Y, B = X 
for A = X,B = Y 



(A10) 



(All) 



spin-orbit 

S 2 E 



53 q ,{v')5 Pq {r) 
S 2 E 
Sj k . q ,(r')Ssi- q (r) 

and tensor terms 
S 2 E 



5J q ,(r')SJ q {r) 

5 2 E 
ST q ,(r>)6s q (r) 



= (& 4 + &iV)Vrtf(r-r / ), (A15) 
= (6 4 + b' 4 5 qq ,)(e k i m V mi r)S{r - r') , 



= -2(6i+SiV)<Kr-r') , (A16) 
= (6!+6iV) 5 ( r - r ') (A-17) 



are most important. 

SRPA equations presented above are obtained for ar- 
bitrary functionals -E[J^(r, i)], including Skyrme ones. 
The model is self-consistent in the sense that both the 
static mean field 



aq y 



(A18) 



and the residual interaction ()A1|) . (|A4|) - (| A7|l are derived 
from the same functional. The rank of the RPA ma- 
trix (|A8p ~ (|A9p is determined by the number K of the in- 
put operators in (|A2[) . As was mentioned above, us ually 
K = 2 -r 5 and so the rank is very small [H, H3, H, [Hf . 
This drastically simplifies RPA computational effort and 
allows to perform systematic explorations even for heavy 
deformed nuclei. 



Here (ph\A q , k ,\) is the matrix element for the two- 
quasiparticle state \ph), e p h is the energy of this state, 
u> v is the energy of the RPA state \v). The amplitudes 
of the RPA phonon operator 



^ ~ E E \ c ph Kh ~ c ph A?h 

q ph£q 

are determined via solutions of (IA8D-(IA9D 



^ph<£q 



(A12) 



? q , k ,lyh\X q q .w) T ifyu (ph\Y q \,) 



q' k' 



2(e p h ± u v ) 



(A13) 



and A h and A p h are operators of creation and destruc- 
tion of two-quasiparticle states. 

Following (|A"4|l - (|A"7|l . the separable ansatz (|A"Tj) ex- 
plores the residual interaction of the Skyrme functional 
through the second functional derivatives. The calcula- 
tions show that for spin-flip magnetic modes the spin 



S 2 E 



8s q >(r')5s q (r) 



fco-fc + ^"(r) (A14) 



J%p°(r)V - (b 2 -b' 2 5 qq ,)A r 



S(r - r') , 



APPENDIX B: SRPA STRENGTH FUNCTION 

Giant resonances in heavy nuclei are formed by many 
RPA states whose detailed contributions cannot be re- 
solved experimentally. Then, instead of solution of Eqs. 
(|A8[) - (IA9[) , a direct computation of the strength function 
(U|) is more efficient and reasonable. In SRPA it reads 



(Bl) 



- z=tj+iA/2 



S{M1-lo) = ^|(^|M|*o>| 2 C(^-^) 

z L Ep^F 0f3 ,(z)D (z)D f3 ,(z) 
ttF(z) 

+ E E £ph(ph\M\) 2 ((u - e ph ) 

q phe q 



where /3 = qkr with r being the time parity, 3 means 
the imaginary part of the value inside the brackets, F(z) 
is the determinant of the RPA matrix (|A"8] ) -(|A"9"] ) with 
lu u replaced by the complex argument z, Fppi{z) is the 
algebraic supplement of the determinant, and 



j>8?(*)=EE 

q' phEq' 



u u (ph\X« k \)(ph\M\) 



4h - z2 



(B2a) 



7 



q' phQq 1 



- ph (ph\Y Q i\)(ph\M\) 



(B2b) 



In the present study we are interested only in spin-flip 
Ml mode. So the strength function is calculated only for 
fx = 1 branch of magnetic dipole excitations. The branch 
with fi = is omitted as it does not support pure spin-flip 
transitions but only its mixture with orbital modes. 



APPENDIX C: INPUT OPERATORS 

The SRPA formalism itself does not provide the input 
operators P q k( r ) in the scaling transformation (|A2|) . At 
the same time, their choice is crucial to converge the 
approximate residual interaction to the true Skyrme 
one with a minimal number of the separable terms. We 
achieve this aim by using P g fc(r)-inputs which compel the 
separable operators X qk {r) and Y qk {r) to have maxima 
in different spatial regions of the nucleus, both in the 
surface and interior. The analysis shows that this way 
indeed allows to get good convergence already with a 
few separable terms. 

The physical arguments suggest that the leading scal- 
ing operator P q i(r) should have the form of the applied 
external field in the long-wave approximation, in our case 



of magnetic field of multipolarity Xfi = 11. In present cal- 
culations we take decoupled spin and orbital input oper- 
ators Pqi = s q + and P q 2 — l\ and supplement them by 
P q 3 = r 2 s q + and P q 4 = r 2 l q + . Then altogether we have 
K=4 input operators and the corresponding operators 
X qk (r) and Y qk (r) have maxima in both surface and in- 
terior of the nucleus. 

In deformed nuclei we should take into account the 
coupling between magnetic and electric = 1 + states. 
So the quadrupole input operator Q q $ — r 2 Y2i(@) with 
the counterpart P q § — i[H, Qgs] is added. Then K=5 
and we have the RPA matrix of the rank 4K=20. 

In terms of two-quasiparticle matrix elements, the re- 
lations between input operators are reduced to 

Q qk (r) = i{H,P qk ] -» 
(ph\Q qk \0) = 2e ph ( P h\P qk \0) - { P h\X q qk \0) (CI) 

for magnetic modes and 

P qk (r) = i[H 7 Q qk ] -> 
(ph\P qk \0) = 2e ph ( P h\Q qk \0) - { P h\Y q k \0} (C2) 

for electric modes. 

For more details see [H, S3, U HI • 



[1] M. Bender, P.-H. Heenen, and P.-G. Reinhard, Rev. Mod. 

Phys. 75, 121 (2003). 
[2] D. Vretenar, A.V. Afanasjev, G.A. Lalazissis, and P. 

Ring, Phys. Rep. 409, 101 (2005). 
[3] J. R. Stone and P.-G. Reinhard, Prog. Part. Nucl. Phys. 

58, 587 (2007). 
[4] P.-G. Reinhard, Nucl. Phys. A649, 305c (1999). 
[5] P. Klupfel, P.-G. Reinhard, T. J. Biirvenich, 

J. A. Maruhn, Phys. Rev. C 79, 034310 (2009). 
[6] T. Lesinski, M. Bender, K. Bennaceur, T. Duguet, and 

J. Meyer, Phys. Rev. C 76, 141312 (2007). 
[7] M. Bender, J. Dobaczewski, J. Engel, and W. 

Nazarewicz, Phys. Rev. C 65, 054322 (2002). 
[8] S. Fracasso and G. Colo, Phys. Rev. C 76, 144307 (2007). 
[9] P. Sarriguren, E. Moya de Guerra, and A. Escuderos, 
Nucl. Phys. A69, 631 (2001). 
[10] Int. Rew. Nucl. Phys. 7, ed. J. Speth (World Sientific, 

Singapore, 1991). 
[11] F. Osterfeld, Rev. Mod. Phys. 64, 491 (1992). 
[12] M.N. Harakeh and A. van der Woude, Giant Resonances 

(Clarendon Press, Oxford, 2001). 
[13] P. Sarriguren, E. Moya de Guerra, and R. Nojarov, Phys. 

Rev. C 54, 690 (1996). 
[14] RR. Hilton, W. Hobenberger, and P. Ring, Eur. Phys. 

J. Al, 257 (1998). 
[15] H.L. Wortche, Ph.D. thesis, Technischen Hochschule 

Darmstadt, Germany, 1994. 
[16] D. Frekers et al, Phys. Lett. B244, 178 (1990). 
[17] S. K. Nanda et al, Phys. Rev. C 29, 660 (1984); A. 
Richter, Phys. Scripta, T5, 63 (1983). 



[18] R. M. Laszewski, R. Alarcon, D. S. Dale, S. D. Hoblit, 

Phys. Rev. Let. 61, 1710 (1988). 
[19] T. Shizuma et al, Phys. Rev. C 78, 061303 (2008). 
[20] F. Tondeur, M. Brack, M. Farine, and J.M. Pearson, 

Nucl. Phys. A420, 297 (1984). 
[21] P.-G. Reinhard, D.J. Dean, W. Nazarewicz, J. 

Dobaczewski, J. A. Maruhn, and MR. Strayer, Phys. 

Rev. C 60, 14316 (1999). 
[22] N. Van Giai and H. Sagawa, Phys. Lett. B 106, 379 

(1981). 

[23] J. Bartel, P. Quentin, M. Brack, C. Guet, and H.-B. 

Haakansson, Nucl. Phys. A386, 79 (1982). 
[24] E. Chabanat, P. Bonche, P. Haensel, J. Meyer, and R. 

Schaeffer, Nucl. Phys. A627, 710 (1997). 
[25] P.-G. Reinhard and H. Flocard, Nucl. Phys. A584, 467 

(1995). 

[26] V.O. Nesterenko, W. Kleinig, J. Kvasil, P. Vesely, P- 
G. Reinhard, and D.S. Dolci, Phys. Rev. C 74, 064306 
(2006). 

[27] V.O. Nesterenko, W. Kleinig, J. Kvasil, P. Vesely, and 

P.-G. Reinhard, Int. J. Mod. Phys. E 16, 624 (2007). 
[28] V.O. Nesterenko, W. Kleinig, J. Kvasil, P. Vesely, and 

P.-G. Reinhard, Int. J. Mod. Phys. E 17, 89, (2008). 
[29] W. Kleinig, V.O. Nesterenko, J. Kvasil, P.-G. Reinhard, 

and P. Vesely, Phys. Rev. C 78, 044313 (2008). 
[30] V.O. Nesterenko, J. Kvasil, and P.-G. Reinhard, Phys. 

Rev. C 66, 044307 (2002). 
[31] V.O. Nesterenko, J. K vasil, W. Kleinig, P.-G. Reinhard, 

and D.S. Dolci, arXiv: |nucl-th/0512045l 
[32] P. Vesely, Ph. D. thesis, Charles University in Prague, 



Czech Rep., 2009. 



